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Abstract 



Since an acceptable dark matter candidate may interact only weakly with ordinary matter 
and radiation, it is of interest to consider the limiting case where the dark matter interacts 
^ ■ only with gravity and itself, the matter originating by the gravitational particle production at 
the end of inflation. We use the bounds on the present dark mass density and the measured 
large-scale fluctuations in the thermal cosmic background radiation to constrain the two 



parameters in a self-interaction potential that is a sum of quadratic and quartic terms in 



a single scalar dark matter field that is minimally coupled to gravity. In quintessential 
inflation, where the temperature at the end of inflation is relatively low, the field starts 



acting like cold dark matter relatively late, shortly before the epoch of equal mass densities 
in matter and radiation. This could have observable consequences for galaxy formation. We 
respond to recent criticisms of the quintessential inflation scenario, since these issues also 
apply to elements of the noninteracting dark matter picture. 



^ ■ I. INTRODUCTION 



We have improving observational evidence [p]] on how structure formed on the scale of 
galaxies and larger from measurements of the angular distributions of the radiation back- 
grounds, surveys of the spatial distributions and motions of the galaxies, and observations of 
the evolution of galaxies and the intergalactic medium back to redshift z ~ 5. It may prove 
useful to complement these advances with explorations of the options conventional physics 
offers for theories of structure formation. Here we consider the possibility that the dark 
matter that is thought to dominate the gravitational growth of structure interacts only with 
itself and gravity, the dark matter originating by gravitational particle production, that is, 
as a squeezed state, at the end of inflation.^ 



Such dark matter particles might be called cabots, in honour of the Cabot family of Boston, who 
were said to be so highly placed in society as to speak only to God. 
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After specifying our model in Section II, we discuss the state of the dark matter field at 
the end of inflation. We begin in Section III with the equilibrium properties of a noninteract- 
ing dark mass scalar field in de Sitter spacetime. For many purposes de Sitter spacetime is 
a useful approximation to inflation, and it has the particular advantage that the equilibrium 
statistical properties of the dark mass field are well specified [@-f§]. However, it is also of 
interest to analyze the effect of the rolling value of the Hubble parameter in a model for the 
relatively late stages of inflation. This is presented in Section IV, following in part the early 
work by Kofman and Linde ||. Section V deals with the constraints on the mass of the dark 
matter particles and on their self-coupling from the conditions that a single scalar field has 
the wanted present mean mass density and field fluctuations that are compatible with the 
large-scale anisotropy of the thermal background. In quintessential inflation |J this dark 
matter candidate has an observationally interesting effective Jeans length. 

We discuss in section VI issues of internal consistency of the model. In a recent preprint 
Felder, Kofman and Linde [0] have presented a stimulating list of potential problems with 
quintessential inflation and indirectly with the noninteracting dark matter model. We ad- 
dress their points in this section. An assessment of consistency with the full suite of obser- 
vational tests requires computations that go well beyond our exploratory discussion; we also 
comment on these open issues in Section VI. 



II. ASSUMPTIONS 

The dark matter is modeled as a single scalar field, y(x,t), that is minimally coupled to 
gravity and has the self-interaction potential 

V = Xy'/i + fiV/2. (1) 

The generalization to several field components with O(N) symmetry and the opposite sign 
of the quadratic part of the potential is trivial and not interesting for the present purpose. 

We assume that there is a time before the end of inflation when the Hubble parameter 
H(t) is close enough to constant that the dark matter field can relax to near statistical 
equilibrium between the amplitude growth driven by quantum fluctuations and the classical 
slow roll due to the self-interaction potential, in an approximation to a de Sitter-invariant 
quantum state. For definiteness, much of our discussion of this assumption is based on a 
quartic inflaton potential, 

U = A^ 4 /4, (2) 

with a self-coupling satisfying ^$ 10~ 14 , for consistency with the isotropy of the thermal 
cosmic background radiation (the CBR). Kofman and Linde H consider the two-field model 
in Eqs (|l|) and (§) in the limit where the evolution of y is essentially classical. Following 
their methods in part, we show in Section VI that, when 

A > A , (3) 

the y-field energy density indeed remains subdominant during inflation, as required for 
internal consistency. We show in Section V that the condition (|3|) is within the ranges of 
values of A and A^ allowed by the CBR. 
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We assume the mass parameter fi is small enough that the quartic part of the potential 
dominates during inflation, that A is small enough that the field may be approximated as 
the sum of a classical part and a free quantum field, and that radiative corrections to the 
potential in Eq. (D) are negligible. The condition for the last assumption is |§ 

3A lni«l. (4) 



64tt 2 A 

We use a homogeneous cosmological model in which the field equation is 

d 2 y ady V 2 y 2 2 

W + 3 am=^- {Xy + » )y - (5) 

For \y 2 <C H 2 , the oscillation frequency due to the quartic part of the potential is small 
compared to H, and we can use the slow roll approximation to the field equation, 

3Hdy/dt = -Ay 3 . (6) 

This neglects the mass term, which is justified for y 2 3> yU 2 /A. 

We consider both the conventional model for reheating and the much lower temperature 
at the end of inflation implied by the quintessential inflation model ||. In the latter case 
the mass density in interacting matter at the end of inflation is proportional to the fourth 
power of the Hubble parameter, with a coefficient that depends on the nature of the matter 
and its interactions. We comment on the uncertainty in this coefficient in Section VI. 

We write the Planck mass as m p \ = [3/(8ttG)}^ 2 in units where h = 1 = c. The Hubble 
parameter during inflation is H and the value at the end of inflation is H x . In section V the 
present value of Hubble's constant is H = 100/i km s" 1 Mpc -1 . Numerical examples assume 
h = 0.7, density parameter Q m = 0.3 in matter capable of clustering, and a cosmologically 
flat universe. 



III. STATISTICS OF THE MATTER FIELD IN DE SITTER SPACETIME 

On scales larger than the de Sitter horizon, H~ 1 , the dynamics of the field y can be 
pictured as a "random walk" [superimposed on the classical slow roll in Eq. (^])] in which y 
undergoes random steps of rms magnitude (5y) TWLS = H/2n per expansion time 5t = H -1 , 
independently in each horizon-size region. The statistical properties of the field y resulting 
from this random process can be described in terms of the probability distribution function 
P(y,t) which satisfies the Fokker-Planck equation 0|| 

9P(v,t) 1 d H 3 d 2 P(y,t) 

= 3Hdy- [V iy)PM] + 8^^^' (?) 

Any initial distribution, P(y, 0), approaches the stationary solution of (0), 

Po(y) = N-'exp(-^-V(y)), (8) 



3H A 



on the timescale 
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r rel ~ \-V*H-\ 



(9) 



The coefficient A" in (§) is determined from the normalization condition f^Po(y)dy = 1. 
The stationary distribution Pq(j/) corresponds to the de Sitter-invariant quantum state of the 
field y. The statistical properties of this state in the case of a quartic potential V(y) = Xy 4 /4 
have been studied by Starobinsky and Yokoyama j|] (hereafter SY). 

The distribution function Po(y) is not Gaussian; in particular it has negative excess 
kurtosis M (</> 4 )/((0 2 ) 2 ) — 3 ~ 0.812. In the following numerical example of correlation 
functions of powers of the field, we will refer to the moments 

(y 2 ) = 0.1318# 2 A- 1/2 , (y 4 ) - (y 2 ) 2 = 0.0206F 4 A" 1 , 

(y 6 ) = 0.01502# 6 A" 3/2 , (y 8 ) - (y 4 ) 2 = .00577 H 8 A" 2 . (10) 
(y 12 ) - 3(y 4 )(y 8 ) + 2(y 4 ) 3 = 0.001755F 12 A~ 3 . 

SY use the Fokker-Planck formalism to compute the field correlation function (2/12/2)- 
Their method and numerical results are readily adapted to get the correlation function of 
powers of the field, in particular the correlation function of the mass density Xy 4 / 4. This is 
done in the Appendix; we find that for a positive integer power n of the field the equal time 
reduced correlation function of y n at large separation ar 12 is 

CnM = fe 2 n > - (y n ) 2 oc x u p , (11) 

where the power law index is 

p = 0.178A 1/2 (odd n), p = 0.579A 1/2 (even n). (12) 

The power law in Eq. (|TTD applies at separations large compared to the comoving coherence 
length 

R c = (aHy 1 exp(l/p) ~ (aH)- 1 e aHTr < i , (13) 

with the appropriate value of p from Eq. (0), and where a is a number of order unity. 
This expression can be interpreted to mean that in the relaxation time r rc i (Eq. |9|) comov- 
ing positions initially separated by the Hubble length have moved to separation aR c , with 
significant loss of memory of the common initial field values at the two positions. 

As an alternative to the Fokker-Planck formalism, it may be convenient to use a discrete 
version of the Langevin equation. In an expansion time, t = H~ l , the value of the field y at 
a fixed comoving position changes by the amount 

y(j) - vU + 1) = vU) - Mi)7(3# 2 ) + hH/(2n). (u) 

The integer j counts successive expansion times H~ l . The Gaussian normal random variables 
tj have zero mean and are statistically independent, so in particular (ijik) = Sj^. The 
second term on the right-hand side of Eq. (|T4T), which represents the classical force due to 
the potential that tends to drive y toward zero, follows from the slow roll approximation 
in Eq. (|6]). The last term represents the frozen quantum fluctuations added to y in the 
expansion time H~ x . 
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In the stationary de Sitter-invariant state, the mean square value of y is independent of 
time. This condition applied to Eq. (0), and keeping first powers of the classical drift part 
and second powers of the fluctuating part, yields the mean mass density in the y-field, 

<P,) = A(?/ 4 >/4 = 3tf 4 /(327r 2 ). (15) 
The condition (y(j) 6 ) = (y(j + l) 6 } similarly yields 

(y 8 )/(y 4 ) 2 = 5. (16) 



These relations agree with the Fokker-Planck results given above, as expected. Eq. (16) 
says the standard deviation of the frozen field mass fluctuations is Sp/p = 2, with coherence 
length R c ~ F- 1 exp(1.7A- 1/2 ). 

Eq. ( |H|) may be applied as a numerical iteration prescription, where the integers j count 
iterations and the numbers ij are drawn from a generator of pseudo- Gaussian independent 
random numbers with zero mean and unit variance. We use it to calculate the two-point 



functions (|TT|) and illustrate the different power law indices for even and odd n. The com- 
putation starts with an initial value y(0) obtained after 300 iterations of Eq. (|T4p starting 
from y = 0. Two time series, yi{j) and y 2 (j), are computed from the same initial value 
y(0) and different sets of the Lj. This is a numerical realization of the values of the y-field 
at successive expansion times H~ l at two fixed comoving positions that are separated by 
the Hubble length at j = 0. A set of realizations is obtained by using the last value of the 
time series y\ as the initial value for the next realization of yi{j) and y 2 (j)- The mean of 
yi(j)y2(j) n across a set of these realizations is an estimate of the expectation value (y^y^) 
at separation xi 2 oc e? . 

Fig. 1 shows numerical results for n = 1, 2, 3, and 4 in Eq. fllTI). The variables are scaled 

to 

c n = X n / 2 H- 2n c n , x 12 = {Hax l2 ) V ~\ (17) 

to remove the dependence on the parameters H and A at zero separation (Eq. pi]) and at 
large separation (Eq. |12| ). The heavy lines in the figure are the averages across realizations 
for two parameter choices, A = 0.1 and A = 0.0001. The latter requires a much larger number 
of iterations to reach a given range of values of the correlation functions. That reduces the 
number of realizations in the computation, so the numerical noise is larger. The straight 
dotted lines are interpolations based on the SY results: the intercepts at zero separation are 
the one-point moments in Eq. (pi]) and the slopes are the power law indices 0.178 and 0.579 
for the scaled variable (Hax)^ (Eq. [12]). Within the fluctuations from the limited number 
of realizations in the averages (1 x 10 8 for A = 0.1 and 1 x 10 7 for A = 1 x 10~ 4 ) the numerical 
estimates are consistent with the expected scaling with A. Since the power-law asymptotics 
apply only for x > R c , one could not expect them to match with the one-point moments 
which correspond to ax ~ H^ 1 <C aR c . Fig. 1 indicates, however, that the difference from 
a pure power law is not dramatic: the mismatch is no more than a factor of about 2. 

The lowest curves in Fig. 1 show the mass correlation function £(£12) oc (yfy 2 ) — (y 4 ) 2 - 
At large separation the numerical result is a factor of about two below the extrapolation 
from zero separation. Since the reduced second moment of the mass distribution is £ = 4 at 
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zero separation (Eq. [16]), a useful approximation to the dimensionless dark mass correlation 
function is 



£(x) = (8 1 5 2 ) ~ 2(x iE /x) £ , e = 0.58A 1 / 2 , (18) 
where x x is the comoving Hubble length at the end of inflation and the density contrast is 

5 = p/( P )-l = yV(y A )-l. (19) 



Comparison to the moments of large-scale fluctuations in the counts of galaxies flO 
requires higher moments of the ^/-matter distribution; we present the example of the three- 
point function. The reduced dimensionless three-point mass correlation function is 

= {vhAvtiltf? - ((ytvt) + (vtvt) + (vivt)) I V> 2 + 2. (20) 

The arguments of £3 are the lengths of the sides of the triangle defined by the positions at 
which the field is evaluated. It is shown in the Appendix that if the three points define 
an equilateral triangle with side x much larger than R c then the three-point function varies 
with x as £3 oc £ 3//2 , where the two-point function £ (Eq. |T3|) is evaluated at x. The one-point 
moments (Eq. pi]) show that £s/£ 3 / 2 = 4 at zero separation. 

Fig. 2 shows numerical realizations of £s/£ 3//2 for equilateral triangles. The estimates of 
£ 3 are based on sets of three independent time series yi(j), 2/2 (j); and y%(j) with the same 
initial value. The values of A, the numbers of realizations, and the scaled variables are the 
same as for Fig. 1. The curves are close to the SY one-point moments at small separation, 
and, as also predicted by the Fokker-Planck method, at scaled separation x ^ 10 the ratio is 
close to independent of triangle size, within the considerable noise from the limited number 
of realizations, at 

c 3 /e /2 - 0.8 (21) 

This gives a useful working approximation to the normalization of the three-point function 
on scales of interest for structure formation.0 

The scaling of the n-point mass correlation functions with distance x at fixed values of the 
ratios x^/xki of distances among the n points may be obtained by the following argument. 
Let 5i = 5p/ p be the mass density contrast at expansion parameter a = ai smoothed within 
a window of fixed comoving size, shape, and position, and let 5f be the density contrast 
smoothed in the same window at a later time, at expansion parameter ay. A realization S(t) 
of the evolution from 5i to 5f is a result of the process of freezing of quantum fluctuations and 
the classical evolution toward 5 = 0. If aj/ai is large enough that the epochs are separated 
by many relaxation times, the initial value 5i is a small perturbation to this process, and 



2 The asymptotic value of the ratio in Eq. (21) can be found by a numerical computation of the 
eigenfunction $2(2/) and the integral C in Eq. ( 110 ) of the Appendix. We have not attempted this 
computation. 
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in the lowest nontrivial order in perturbation theory the expectation value of 5f for given 
initial value Si is linear in 8f 

(Sfh = (22) 

The transfer coefficient T is a function of the window size and shape, which we are holding 
fixed in comoving coordinates, and of the expansion factor a//aj, which we will vary. Next we 
consider two disjoint windows, both of fixed size, shape and position in the same comoving 
coordinate system. If the distance between the windows is much larger than the Hubble 
length at the evolution of the density contrast in window 1 is independent of the process 
in window 2, and it follows that the equal time mass two-point correlation function satisfies 

(*,(1)*,(2)> =r a G(l)tf 4 (2)>. (23) 

In the de Sitter equilibrium state the two-point mass correlation function depends only on 
the proper separation ax, so we can rewrite Eq. ( p3| ) as 

£(6z) = T{bfi(x) } b = af/ai. (24) 

This implies 

f oc x~ 2u , T oc b~ v , (25) 

where v is a constant that must be positive for convergence. That is, at large separation x 
the correlation function varies as a power of x, as was shown more directly by our use of the 
SY analysis. The n-point function similarly satisfies 

(8 f (l)...8 f (n))=T n (S i (l)...S i (n)), (26) 

from which it follows that the n-point function scales with the separation x at fixed ratios 
of separations among the n points as 

Cn oc x- w , (27) 

consistent with Eq. (pip . 

Eq. (p7|) says the n th central moment of the dark mass M contained in a volume of 
fixed shape scales as M n = ((M - (M)) n ) oc a n , where a = (((M - (M)) 2 )) 1 / 2 is the 
standard deviation (and assuming the integrals over the correlation functions converge at 
vanishing separation between points). The smoothed mass fluctuations do not approach a 
Gaussian as the size of the smoothing window is increased. Rather, this is a scale-invariant 
fractal in which the mass fluctuations smoothed and referred to the standard deviation are 
independent of the size of the smoothing window. 



IV. THE DARK MATTER AT THE END OF INFLATION 

We have assumed so far that the evolution of the inflaton <fr is sufficiently slow that the 
field y has enough time to settle into its de Sitter-invariant "equilibrium" state. Here we 
specify the conditions of applicability of this picture under a specific inflaton potential, and 
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we consider the transition to the opposite regime, which applies toward the end of inflation, 
when the evolution of y is almost classical. 

For the inflaton potential in Eq. (0) an approximate solution for the expansion history 
through inflation is Jll| 

H = H x e- H *\ H x = ^\l /2 m ph <\> x = (8/3) 1/2 m ph (28) 



with expansion parameter 



loga x /a = H/H x -1. (29) 



We write the comoving length scale x of field fluctuations frozen at time t, when the 
expansion and Hubble parameters are a and H, as 

Hax = 1. (30) 



This comoving length reaches its minimum value, x x , at the end of inflation. In Eq. ( |28|) 
inflation ends at t = 0, when the Hubble parameter and inflaton field values are H x and 4> x . 

In the approximation of Eqs. ( f28| ) to ( |30| ) the Hubble parameter H when the field 
fluctuations on the comoving length scale x are freezing satisfies 

\og(x/x x ) ~ hg(H x /H) + H/H x - 1. (31) 

A useful approximation at a C a x is 

x/x x ~e H/H *. (32) 
The relaxation time for y and the evolution time for H are (Eqs. || and |28| ) 

T re l ~ \~ l/2 H~ l , T H ~H-\ (33) 

When r ro i th the y-field relaxes to statistical equilibrium between fluctuations and clas- 
sical slow roll, in a good approximation to the de Sitter equilibrium state. When r re \ ^> th 
the field evolution is close to classical. At the transition between these limiting cases, at 
expansion parameter a e , the Hubble parameter is 

H e ~ \~ 1,2 H X . (34) 

Since the y-field at a = a e is close to statistical equilibrium its characteristic value is 
y e ~ A~ 1//4 if e (Eq. |1^). This is related to the value of the inflaton at a = a e by 

yj^ ~ (A^/A) 1 / 2 . (35) 

Kofman and Linde [BJ show that, if the inflaton and the y- field both have quartic potentials, 
the evolution in the slow roll classical approximation is 

2\-l f\„.2\ 



- 1 - (Xy 2 y l = constant = (A^)~ - (\y x Y . (36) 

Eq. (|3~5|) says the constant of integration is not very important at a = a e , so a good 
approximation to the smaller field values at a e C a ^ a x is 
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y(t) = (X^/Xy^it). 



(37) 



We consider first the field fluctuations frozen at a e < a < a x , when the evolution of y is 
close to classical. In unit logarithmic interval of x the contribution to the variance of the 
frozen dark mass field is 5y% = (H/2ir) 2 . Apart from regions where y happens to be close to 
zero, this is a small fractional perturbation to the field value. With Eqs. fl3"6| ) and Q3"7| ) we 
see that the field perturbation at the end of inflation isQ 



Vx 




(38) 



Here H, y, and <ft are evaluated when field fluctuations on the scale x are frozen, at 
log(x/x x ) = (4>/4> x ) 2 in the approximation of Eq. (|32|) . With Eq. (p8|) we get 

Sy x /y x ~ 0.1A 1 / 2 [log(x/x !r )]- 1 / 2 . (39) 

Since this is nearly independent of the length scale x the power spectrum varies with wave 
number about as k~ 3 , and the mass correlation function is 

£~ MAlog(x e /x), (40) 

where M is a dimensionless constant. 

At separations larger than the field coherence length, R c ~ (a e H e ) ~ 1 exp(l/p) with p ~ 
A 1 / 2 (Eq. |T^), Eq. (fh]) fails; the values y e may be quite different at the two points. The 
classical evolution Eq. (|36|) says the field values at a given comoving position a = a e and at 
the end of inflation at a = a x satisfy 

iH^ 1 -^- (41) 

If A <C 1 this says the field value is y x ~ ±X~ l l 2 H x everywhere except near the surfaces 
where y e vanishes, and the energy density is 

Pm - ^~ l H x , (42) 

except near the zeros of y e . If the distance between zeros were larger than the typical motion 
of the dark mass after inflation the present distribution would be close to uniform apart from 
surfaces of low density. 

We estimate the form of the mass correlation function at the end of inflation and for 
separation x 3> R c by the argument used to obtain the scaling of the n-point mass correlation 
functions in the de Sitter-invariant case in Eq. ([27|). Consider two windows with fixed sizes 
and positions in comoving coordinates. The dark mass density contrasts smoothed within 



3 To avoid confusion we remind the reader that Sy x /y x is the fractional fluctuation in the dark 
matter field value at the end of inflation. The subscript x in Eq. ( |38| ) has nothing to do with the 
comoving length scale of the fluctuation in Eq. (p9|). 
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these windows are 4(1) and 4(2) at expansion parameter a = a e , when the field fluctuations 
start to depart from the de Sitter-invariant state, and 4(1) and 4(2) at the end of inflation 
at a = a x . If the comoving separation a e X\2 of the windows is large compared to R c at a e 
then we have, following the derivation of Eq. ( p5|) . 

(4(1)42)) = T 2 (4(l)4(2)) ~ N{x x /x 12 f 6VX ~ Cx{x 12 ). (43) 

The power law coefficient follows from the de Sitter- invariant correlation function ([H|). The 
transfer coefficient T, defined as in Eq. fl2"2p, depends on A. In general it also depends on 
the window size and shape, but in the limit where the two windows are small compared to 
their separation Eq fl43|) reduces to the two-point mass correlation function. 
Quite similarly, the three-point mass correlation function satisfies 

(5(1)4(2)4(3)) = T 3 (4(l)4(2)4(3)). (44) 



Combined with Eq. ( |43| ) this implies that the relation ( pT|) between the three- and two-point 
correlation functions is still valid at the end of inflation. 

If A is not much smaller than unity the near de Sitter evolution ends not long before 
the end of inflation. In this case an estimate of the mass correlation function at the end 
of inflation from a numerical realization of the process is easy and useful. As in Section II, 
we label successive e-foldings of the comoving Hubble length by the integer j, but now j 
decreases with time, with j = 1 at the end of inflation. The values of the Hubble constant 
follow from Eq fl31"|), 

r- pHj/H x -i 

J, W («> 

Eq. (|36| ) says the values of the dark matter and inflation fields at successive e-foldings are 
related by 



-(--- 

A \y] Vj+iJ h \(f>j "j + i . 

2 




(46) 



The result of multiplying this by H£, rearranging, and adding the quantum noise term 
±H/ (27r) at each e-folding of x is 

where the Lj are independent Gaussian normal numbers, the field has been scaled to 

z = y/H x , (48) 

and 
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The Hj come from the numerical solution of Eq. (|45|). The two-point function at separation 
x/x x = e- 7 is the mean of products of field values computed starting from a common value 
j time steps before the end of inflation. 

Figure 3 compares the mass correlation functions at the end of inflation for constant 
Hubble parameter and for the rolling case in Eq. fl28|), for A = 0.1. One sees that at 
relatively small separations £(x) is considerably flatter in the rolling H expected 
from Eq (f|0p. The values and rates of change of the correlation functions are roughly similar 
at separations x ~ 10 15 x x .. This suggests that for A = 0.1 the constant N in Eq. fl43|) is not 
greatly different from unity. 

At smaller A a numerical realization of the two-point function is difficult because the 
relaxation time is long, but it is easy to get the one-point distribution of y at the end of 
inflation. At A = 0.001 the distribution is strongly peaked at y ~ iA -1 / 2 ^, as expected 
(Eq. |42| ), but there is significant scatter from the quantum fluctuations appearing at a ^ a e . 
For A = 0.1 the distribution of field values y is peaked at zero, and not greatly different from 
the stationary de Sitter case. 

In the next section we use the estimates of the field mass fluctuations in Eqs. (p9|) and ( fj3|) 
with the measured CBR anisotropy to find bounds on the parameter A. We will argue that 
the considerable uncertainty in N translates to a relatively small uncertainty in the bound 
on A. 

V. THE DARK MATTER AT THE PRESENT EPOCH 

We consider three constraints: the present mean mass density in dark matter, the small- 
scale cutoff in the gravitational growth of clustering of the y-matter, and the large-scale 
anisotropy of the thermal background radiation. Because most of our estimates are quite 
approximate we ignore most numerical factors that are of order unity. 

A. The Mean Mass Density 

With the change from proper world time t to conformal time x° = J dt/a, and the 
definition y = ay, the action for the ?/-field is 

S = Jd'x (lytf* - \\y* - ^<??R) , (50) 

where the index is raised by the Minkowski metric tensor and we ignore the quadratic part of 
the potential. The Ricci tensor R is on the order of the square of the Hubble parameter a/ a 
when the universe is matter-dominated, and it is well below that when radiation-dominated. 
When the Ricci term may be neglected a solution for y in Minkowski spacetime is a solution 
for ay expressed as a function of comoving coordinates and conformal time. The j/-field 
energy density depends on the proper time derivative y = (dy/dx° — ya)/a 2 . If the field 
oscillates rapidly on the scale of the Hubble time, either because the potential is driving 
rapid oscillations or the field varies with position on scales small compared to the Hubble 
length, then y = a~ 2 dy/dx° to good accuracy, and the proper field energy density is well 
approximated as the energy density of y in Minkowski coordinates divided by a 4 . Since the 
former conserves energy the mean energy density in the ?/-field varies as 
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Py oc a 4 . (51) 
The characteristic amplitude of the field oscillations thus scales as 

y<xl/a(t). (52) 



Ford [12] obtained these results for a quartic potential when y is a function only of world 
time. It is an easy exercise to obtain Eqs ([51]) and fl5"2|) by extending Ford's method to the 
case where the scales of length and time variations of y both are small compared to H" 1 . 
The mean field value at the end of inflation is, from Eqs. (^) and Q3"7|), 

y x ~ \~ 1/2 H XJ (53) 

At this field value the characteristic field oscillation frequency is lo ~ A 1 / 2 ^ ~ H x . This 
means the field starts oscillating at about the time inflation ends, so the field amplitude after 
inflation decays as y ~ y x a x /a until the quadratic part of the potential becomes comparable 
to the quartic part, at y ~ y^ = /iA -1 / 2 . The expansion factor at which this happens is 

a^/a x ~ H x //i. (54) 

The present mean mass density in the y-field, which we are assuming is comparable to the 
total, is 

p ~ z eq T^ ~ A~ 1 iJ^(a x /a At ) 4 (a M /a ) 3 . (55) 

The mass density at the end of inflation is A?/ 4 /4 ~ A _1 if 4 (Eq |53|). The present CBR 
temperature is T Q , and the redshift at equality of mass densities in matter and radiation is 
z eq ~ 3500 (for Vt m = 0.3 and h = 0.7). 

In the conventional model for reheating it will be sufficient to consider the case where 
the expansion is radiation-dominated at temperature T x ~ (m p iif ) 1//2 at the end of inflation. 
Then the expansion factor from the end of inflation to the present is a Q /a x ~ T X /T Q . With 



Eqs. (|5J) and (|55l) we get 

At ~ Xz eq T m 3 pl /T^ ~ 10 5 AT U 3 GeV, (56) 

for temperature T x = 10 14 T 14 GeV at the end of inflation. The redshift at which the dark 
matter starts acting as a massive field is 

aja^ ~ \z cq (m pl /T x )* ~ lO^ATf/. (57) 

If A ^> 10 14 and T 14 ~ 1 the length scales of interest for extragalactic astronomy appear at 
the Hubble length at redshifts well below this value of a Q /a M , and the dynamical behavior 
of the y-matter is not significantly different from the usual cold dark matter. 

In the quintessential picture the mass densities at the end of inflation in gravitationally 
produced interacting matter and noninteracting matter are 

^(matter) = C# 4 , p^dm) ~ A^tf 4 . (58) 

The constant C depends on the matter interactions, as discussed in Section VI-A. Consis- 
tency with the standard model for the origin of the light elements requires that the ratio of 
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these mass densities, / ~ AC, be greater than about / ~ 15 at the time of light element 
production. If this condition is violated the model is excluded. If the model is viable we 
can normalize the parameter constraints to the value of /. The redshift at which the y-field 
starts acting like a nonrelativistic massive field is 

= a Q /a^ ~ fz eq . (59) 

The radiation temperature at the end of inflation is T x ~ C X ^H X ^ so the expansion factor 
to the present is a a /a x ~ C x ^YL x jT . With Eq. Q51D we have 

\i ~ AC 3/ S eq T ~ 1AC 3/4 eV. (60) 



B. The Effective Jeans Length 

In quintessential inflation y starts acting like a massive field at relatively low redshift 
(Eq. |5Sp, and we have to check the effect on the gravitational growth of clustering of the dark 
matter. The y-field fluctuations that appear at the Hubble length well after inflation have 
wavelengths that are much larger than the period of oscillation of y, and the gradient energy 
density (Vy) 2 /2 thus is strongly subdominant to the kinetic and potential energy densities. 
We assume dynamical evolution leaves the gradient energy subdominant. We proceed by 
deriving a virial theorem valid when the field oscillation frequency is much larger than the 



Hubble parameter, following Ford ||12| , and use it to estimate the effective pressure from the 
trace of the stress-energy tensor. 

Neglecting the cosmological expansion, the field equation is 

y- V 2 y + \y 3 + p 2 y = 0. (61) 

The result of multiplying this equation by y, averaging over space, and averaging over a 
time interval that is much longer than the field oscillation time and much shorter than the 
cosmological expansion time is 

(y 2 ) = ((Vy) 2 ) + \(y 4 )+p 2 (y 2 ). (62) 

With this relation the mean energy density is 

<p> = ((Vy) 2 > + 3A(y 4 >/4 + / iV>, (63) 
and the mean of the trace of the stress-energy tensor is 

(T) = (p}-3(p}=v 2 (y 2 ), (64) 
where p is the effective pressure. If the gradient energy is subdominant we have 

<p> i My 4 ) 



(p) 3\(y*)+4p 2 (y 2 )/3- 



(65) 



Properties of y(t) considered as a single oscillator with the equation of motion y = —Xy 3 
are discussed by Greene, Kofman, Linde, and Starobinsky fll3 |. The energy equation for the 
oscillator is 
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y 2 /2 + Ay 4 /4 = Ay 4 /4, 



(66) 



where the amplitude is y , and we can use this to compute time averages of moments of y. 
In particular, 

(y 2 ) 2 /(y 4 ) = 48[r(3/4)/r(l/4)] 4 ~ 0.63. (67) 

When the amplitude is small enough that the quadratic part of the potential dominates, 

(y 2 ) 2 /(y 4 ) = 2/3. (68) 

In short, we can take it that {y 4 ) ~ (y 2 ) 2 , and we see from Eq. (|65| ) that (p) / (p) ~ 1/3 until 
(y 2 ) ~ A* 2 /A, at expansion parameter a ~ a M . In quintessential inflation this is at redshift 
z^ ~ fz cq (Eq. |59|). At z < z^ the ratio varies as (p) / (p) oc a -3 , because (y 2 ) oc a -3 for a 
massive nonrelativistic field. Thus the effective velocity of sound is c s oc a~ 3//2 oc r 3 / 4 at 
z^ > a > z eq and c s oc t _1 at z < z eq . The physical Jeans length ~ c s t increases as t 1//4 at 
z^ z £ z eq and then approaches a constant. The comoving Jeans length c s t/a is maximum 
at Zp, and the maximum comoving Jeans length referred to the present epoch is 

Lj ~ ~ 6/(/fi m /i 2 ) Mpc. (69) 

The world time at redshift z M is the present matter density parameter is Q m , and the 
present Hubble parameter is h. On comoving scales larger than Lj it is a good approximation 
to neglect the effect of the field stress on the gravitational evolution of the y-field mass 
distribution. On smaller scales the dynamical behavior of the field requires a more detailed 
analysis than is attempted here. 

The effective Jeans length in Eq. (|69p is smaller than that of a family of neutrinos with 
mass ~ 30 eV, and may be comparable to the mean distance between large galaxies. 



C. The Large- Scale Thermal Background Anisotropy 

The parameter A is bounded by the effect of isocurvature fluctuations in the |/-field at 
the end of inflation on the angular distribution of the thermal background radiation (the 
CBR). For definiteness we present numerical results for quintessential inflation. 

This analysis of the temperature fluctuations on large angular scales uses the simplifica- 
tion that the radiation pressure gradient force has little effect on the mass distribution on 
the scale of the present Hubble length, so we can imagine that well-separated regions we see 
at the Hubble length evolve as separate homogeneous cosmological models. In one of these 
homogeneous models the field value and radiation temperature at the end of inflation are 
y x and T x . By repeating the computation in Section IV-A one finds that the present dark 
mass density varies with these parameters and the present temperature, T D , as 

Po ~ XWC^iifaTo/T,)*. (70) 

Under the isocurvature initial conditions (and assuming the ?/-mass density at high redshift 
is subdominant, as required for light element production), the temperature T x is nearly 
homogeneous, the same in all the model universes that represent the evolution of different 
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regions. Also, the present dominant mass density, p a , is close to homogeneous on the scale 
of the present Hubble length, because pressure gradient forces cannot generate large-scale 
curvature fluctuations. Thus the present large-scale CBR anisotropy is 

( 7 i) 

T y x 4 Px ' K J 

This relation also follows in the conventional model for reheating. 

We need expansion factors during and after inflation. In quintessential inflation the ratio 
of the co moving Hubble length now, x Q = (HoCio)" 1 , and the Hubble length x x at the end of 
inflation is 

x H x a x sy-i/4. To ,„ n \ 
~ ~ TFT" ~ 6 JF- I 72 ) 

The ratio of the present temperature and Hubble parameter is 

T /H ~ e 67 . (73) 

When field fluctuations during inflation are being frozen at the comoving scale x D of the 
present Hubble length the expansion and Hubble parameters a p and H p satisfy 

H p a p ~ H a ~ C l/A H x a x H /T . (74) 
These relations with Eq. fl29|) say the expansion factor from a p to the end of inflation satisfies 

a x /a p ~ C- l/ \T /H ) \oga x /a p . (75) 

If C is on the order of unity the numerical solution is 

H p /H x = loga x /a p ~ 72. (76) 

We have to consider two cases, where the dark matter field departs from near statistical 
equilibrium and commences classical slow roll before or after ax Q appears at the Hubble 
length during inflation. The Hubble parameter at the breaking of equilibrium is H e ~ 
X~y 2 H x (Eq. |34j) , so a critical parameter value is 

X p = (H x /H p ) 2 ~ 1 x lO" 4 , (77) 

from Eq (|76|) . If A X p then we observe in the CBR anisotropy the effect of fluctuations 
frozen while y was behaving as a classical field that is slightly perturbed by quantum fluctua- 
tions. In this case Eqs. (|39|), (|7lf ) to (f73f) and f[76|) , with the measurement ST /T ~ lx 10~ 5 , 
imply 

A < 1 x KT 6 . (78) 

This assumes C is on the order of unity in Eq. ([72|) , and unless C is very large its value 
does not significantly affect the bound on A. The condition C ~ f/\ from Eq. (|5"S|) with 
/ ~ 10 requires C <L 10 7 . This limit is close enough to unity for the purpose of the numerical 
estimate in Eq. (|7"6f) . 
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The bound in Eq. ( [T8|) differs from the analysis of Felder, Kofman and Linde by the 



factor in Eq. (38) that takes account of the classical decay of y from a p to the end of inflation. 

In the other limiting case, A 3> X p , it is convenient to use the expansion of the observed 
angular distribution of the background temperature in spherical harmonics, with coefficients 

a?= f (KlY^STo/To. (79) 




With Eq (|71]), and using the addition theorem, J2 m Y l m {l)Yf m (2) = {21 + l)P,(cos0i 2 )/47r, 
we can write the mean square value of a multipole moment as 

= ^^^KT> = J ^(cos0)e(*)/16. (80) 

In this normalization, (5Ti/T) 2 is the contribution to the variance of the sky temperature 
per logarithmic interval of the spherical harmonic index /. The argument of the dark mass 
correlation function at the end of inflation is x = 2x sin#/2, where the present angular size 
distance x Q back to high redshift is 

x ~ 3.2CEU,)- 1 . (81) 

Here fl m = 0.3, we assume zero space curvature, and the present Hubble and expansion 
parameters are H Q and a . 

For the dark mass correlation function in Eq. (f43|) the large-scale anisotropy spectrum is 

(y) 2 = ^^y^ © ow 7-i - 2 ^°- 3V ~ x ~ 1 x io "°- ( §2 ) 

The observational bound applies at I ^ 30. We have written x x /x = (x x /x )(x /x); the 
second factor produces the last factor in the integrand. The solution to this equation at 
A <C 1 is not relevant because the limit in Eq. (|78|) applies. In the solution at A close to unity 
the integral is of order unity and the value of the integral and the factor iV in the primeval 
mass correlation function do not much matter because the sensitive function is (x x /x o ) ' 6X . 
We have 

A > 0.3. (83) 

In the conventional reheating picture log a x /a p is about ten percent smaller. To the 
accuracy of our estimates the constraints on A are the same. 



VI. DISCUSSION 

We comment first on issues of consistency and reasonableness of the noninteracting dark 
matter picture within conventional and quintessential endings of inflation, and then take 
note of some observational challenges. 
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A. Theoretical Issues 



For the purpose of exploring a simple example of dark matter that interacts only with 
itself and gravity we have considered a quartic plus quadratic self-interaction potential. This 
functional form has the arguably attractive feature that one can choose constant coefficients 
that imply an interesting present mean dark mass density and CBR anisotropy. We are 
not competent to say whether other models for a self-interaction potential would be more 
plausible from the point of view of fundamental physics or would produce dark matter 
candidates with more interesting properties. 

In the quartic plus quadratic potential the isotropy of the CBR requires that the dimen- 
sionless parameter A be in the narrow range 

0.3 < A < 1, (84) 

or else be quite small, 

KT 6 > A > A < 1(T 14 . (85) 

The last bound is the condition on the adiabatic density fluctuations produced by the infla- 
ton. The upper limit in Eq. (p4j) is the condition that the potential be safe from significant 
renormalization (Eq. |J). We know of no reason to think either range of values of A is partic- 
ularly attractive within fundamental physics, though it is to be hoped that input from this 
direction eventually will be a factor in the completion of a satisfactory theory of structure 
formation. 

In a recent paper Felder, Kofman and Linde JF| (hereafter FKL) raise a number of 
issues relevant to our analysis. They note that scalar fields whose particles are produced 
gravitationally at the end of inflation could end up dominating the energy density, thus 
assuming the role of the inflaton. This could indeed happen if the self-couplings of some 
of these fields in quartic potential models were smaller than that of the inflaton. One sees 
from Eqs. (|35| ) and (|3?D that the ratio of energy densities in the inflaton and a field with 
self-coupling Xj in the quartic model is 

Pj/p<t> = A /A i . (86) 

Thus consistency within this model requires Xj ^> A^. Since the interacting scalar fields of 
the quintessential model are usual spin-0 particles, like the Higgses that interact with 
gauge fields, their couplings are not expected to be particularly small. The inflaton coupling, 
on the other hand, is required to be small [Eq.fl8"o])], and there does not seem to be a problem 
in meeting the conditions on Xj. 

FKL note that the Xj are also constrained by the condition that the isocurvature fluctu- 
ations produced by the scalar fields not violate the CBR isotropy. The point is valid when 
some of the fields represent stable non-interacting massive particles. Indeed, we find that 
the range 10~ 6 <C Xj £ 0.3 is excluded for noninteracting dark matter.^ 



4 Our larger lower bound on \j seems to be mainly due to the correction for the decay of field 
fluctuations in Eq (^). Our larger allowed range applies when the equilibrium between quantum 
excitations and classical decay of the dark matter field persists to close to the end of inflation, a 
case FKL do not consider. 
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Interacting scalar fields, such as the Higgses, also are produced with inhomogeneous 
distributions at the end of inflation. These particles interact by the usual decay, annihi- 
lation and production processes that produce local thermal equilibrium. In the absence of 
particle number conservation laws, this local thermal equilibrium is fully described by a 
single function of position — the temperature or the mass density. Apart from the adiabatic 
perturbations from the inflaton the model predicts negligible spacetime curvature fluctu- 
ations at the end of inflation. This means the matter density fluctuations correspond to 
fluctuations in the local starting times of cosmological expansion with universal values of 
the cosmological parameters. If all matter interacts and relaxes to local thermal equilibrium 
the result at the present time is an unacceptably homogeneous universe. If a dark matter 
component does not interact with the rest of the matter, the possibility considered in this 
paper, fluctuations in composition, which is to say isocurvature fluctuations, remain. There 
may also be adiabatic curvature fluctuations from the inflaton, of course. Either or both 
could act as seeds for structure formation. 

FKL also point out that the moduli problem is more severe than in the conventional 
picture for reheating. Moduli are light scalar fields with only gravitational-strength coupling 
to ordinary matter. They necessarily arise in supergravity models and superstring theories. 
The curvature correction generally makes the effective potential of a modulus \ during 
inflation different from at late times, the minima of the two potentials tending to be displaced 
by Ax ~ m v \. Thus, after inflation one has a nearly homogeneous field % with p x ~ 
m^(Ax) 2 ~ m^mpj. The mass m x is determined by the supersymmetry breaking scale 



Vsusy, 



m 



x 



Vsusy / m P i- ( 87 ) 



In models with gravity-mediated supersymmetry breaking, tjsusy ~ 10 11 GeV and m x ~ 
10 3 GeV. The main problem with the moduli is that they may decay very late. With 



m x ~ 10 GeV, the lifetime is 



m\ x jm\ ~ 10 5 s, (88) 



so one runs into problems with nucleosynthesis and with photodissociation and photopro- 
duction of light elements by the decay products of % [HI • 

There are two commonly discussed possible solutions to the moduli problem. First, a 
short period of secondary inflation could dilute the moduli ||15|| . This naturally occurs in 



"thermal" inflation [[HJ, but seems unlikely to be effective in quintessential inflation, where 
the relative density of the moduli is much higher than in the conventional picture, as pointed 
out by FKL. Second, the minima of V(x) during and after inflation may coincide due to 
some symmetry [17-IIJ. In this case moduli are produced only gravitationally, like other 



light scalars, and have density p x ~ H 4 at the end of inflation. The density of matter in 
standard inflation is p m ~ m^H 2 , and p x j p m ~ H 2 /mh is small enough to solve the moduli 
problem. However, in quintessential inflation p m ~ H and p x j ' p m ~ 1, which is clearly too 
high. 

There are other ways out, however. During inflation, moduli typically acquire masses 
fT7| , |2m m x ~ (3H with /3 > 1, and it is conceivable that /3 > 1 jJT7| , pT| . The gravitational 
production of moduli would then be exponentially suppressed. Moduli can also develop 
non-perturbative potentials, independent of supersymmetry breaking p2 . Their mass can 
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then be large enough that the lifetime in Eq. fl88|) is r Is. Assuming coincident minima of 
V(x) during and after inflation, this would eliminate the problem with light elements. If the 
moduli decay prior to baryogenesis, then the relaxation to thermal equilibrium eliminates the 
isocurvature perturbation associated with the inhomogeneous primeval local ratio of entropy 
density to the number density of modulus quanta. Otherwise, the isocurvature perturbation 
survives in the form of an inhomogeneous ratio of the baryon and photon number densities 
and can later play a role in structure formation. 

Another possibility is that the supersymmetry breaking scale could be much smaller than 
10 11 GeV; in some models it can be as low as the electroweak scale 12311. The moduli would 



then be very light, their lifetime would be much greater than the present age of the universe, 
and problems with nucleosynthesis would not arise. (Here we still assume coincident minima 
and also that the moduli develop masses ~ H during inflation). 

There is a related issue for gravitinos. The lower energy density at the end of inflation in 
the quintessential picture reduces thermal gravitino production but increases the problem 
of gravitational production, as FKL note. Apart from potential energy terms the gravitino 
obeys a conformally invariant field equation |[24|| . If during inflation the gravitino had a 
substantial effective mass, the gravitino energy density from gravitational particle production 
could be unacceptable in quintessential inflation. We are not aware that the effective mass 
has to be this large, however. 

The moduli and gravitino quanta are close approximations to the proposed noninteract- 
ing dark matter, but with potentially unacceptable production of mass density. This is a 
challenge, but within the not inconsiderable uncertainties of a supersymmetric theory that 
has not yet been fully specified the challenge does not seem serious enough to discourage 
further exploration of the picture. 

We note finally that the density of matter produced in quintessential inflation can be 
substantially higher than our original estimate in Ref. ||. The quanta of spin-0 fields Xj 
produced due to the change in the expansion law at the end of inflation have wavelength 
~ H~ x and energy density pj ~ In addition to these short-wavelength quanta, there 



is also a nearly homogeneous component of Xj- The classical solution (|37f ) always applies at 
the end of inflation, and the energy of the homogeneous component follows as in Eq. (|58|), 



>f ~ (Aj/A^mJ, ~ XfHi, (89) 



Pj 

where H x ~ X]^ 2 m p i is the expansion rate at the end of inflation. We see that pf ] > pf for 
small Xj. For ordinary (interacting) matter fields, the homogeneous component decays and 
thermalizes, just like the short-wavelength quanta. We note also that this enhancement of 
the energy density in the homogeneous component disappears if the fields Xj develop masses 
rrij ~ H during inflation. In this case the original estimate of ||, pj ~ is still valid, for 
Xj < 1, because the quartic term is small compared to the quadratic part of the potential. 



B. Observational Issues 



Previous examples of noninteracting dark matter [25], [26|| postulate a quadratic self- 
interaction potential with a time-variable mass. These are CDM models, where the CDM 
is defined as matter that has had negligible pressure since appearance at the Hubble length 
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of all length scales of astrophysical interest. The two-field CDM model in [26[] could be 
adjusted to fit the spectrum STi/T of angular fluctuations of the thermal background, to 
the accuracy of the measurements then available, but at Z ~ 50 the model power is a factor 



of two above the new Python |27j result. Bridging this gap would require even greater 
contrivance. This isocurvature CDM model thus does not seem to be viable, but it does offer 
a useful reference for a first assessment of observational tests of the quartic plus quadratic 
self-interaction model. 

In the isocurvature CDM model the primeval distribution of the radiation is much 
smoother than that of the CDM because the mean mass density in the CDM is much smaller. 
With the quartic potential considered here the dark matter approximates a fluid with the 
equation of state of the radiation until the redshift approaches z eq (Eq. |59|). This means 
a given primeval isocurvature fluctuation Sp y /p y in the dark matter distribution produces 
a much larger primeval perturbation to the distribution of the radiation than for CDM. 
The effect on the fluctuation spectrum of the CBR on angular scales below the limit of 
application of Eq. (^) remains to be computed. 

In the isocurvature CDM model the primeval mass distribution is p oc y(x) 2 , with y a 
random Gaussian process, and the non-Gaussian mass density fluctuations in this model 
may violate the measured skewness and kurtosis of the large-scale fluctuations of galaxy 
counts |T(| . The numerical realizations in Fig. 2 indicate that in the quartic potential model 
the skewness of the primeval mass distribution is a factor of about 3 smaller than in the 
case p oc y 2 . Again, the effect on the observational constraint may be worth investigating. 

In quintessential inflation the dynamical behavior of the y-field would be an impor- 
tant factor in the formation of the first generations of structure (Eq. |69|). Khlebinkov and 



Tkachev |28| present numerical simulations of the three-dimensional behavior of a scalar 
field with quartic self-interaction potential, in connection with the dynamics of reheating in 
conventional inflation. They find a marked tendency for the growth of field fluctuations on 
ever smaller scales. The qualitative effect is easily understood: where the field value \y\ is 
larger than average the field oscillation frequency is larger, and the frequency differences at 
different positions produce growing field gradients. Our preliminary numerical experiments 
in one space dimension plus time suggest the growing field gradients produce energy fluxes 
that redistribute energy and the space distribution of field oscillation frequencies. This tends 
to unwind the field gradients through substantial parts of space, while leaving localized re- 
gions with relatively large field gradient energy density. The effect is curious enough that 
further study would be at least intellectually interesting, and it is conceivable that it adds 
an interesting complication to the observed growth of structure on the scale of galaxies. 

Finally, we return to our opening remark, that it may be useful to complement the 
growing observational basis for a theory of structure formation with surveys of the options 
available within conventional physics. It is natural to consider simplest possibilities first, 
but prudent to bear in mind the possibility of complications. An example is the evidence 
that the density parameter in matter capable of clustering on scales much smaller than the 
Hubble length is significantly below unity If this were valid it would mean that the 
large-scale nature of the universe is not as simple as we could imagine within conventional 
physical ideas, a result that could hardly be described as surprising in light of the ample 
opportunities for complexity in the physical universe. We arrived at a simple model for 
structure formation, adiabatic CDM, fifteen years ago, after abandoning a few even simpler 
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ideas. In the present paper we have considered a class of dark matter candidates whose 
properties could be considerably more complicated than cold dark matter, and the obser- 
vational predictions much harder to work out. On the other hand, the picture does not 
seem unnatural: possibly some of the fields gravitationally produced by inflation are unable 
to thermalize with ordinary matter; perhaps such fields end up with an observationally in- 
teresting dark mass density; perhaps their irregular primeval distributions are a significant 
factor in seeding the gravitational growth of structure. 
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VIII. APPENDIX A. CALCULATION OF TWO- AND THREE-POINT 

FUNCTIONS 

To calculate the two- and three-point correlation functions in de Sitter space, it will be 
convenient to introduce the conditional probability H(y, t\yo, to) for the field to have value 
y at time t, given that it had value yo at time to at the same co-moving position. This 
probability satisfies the Fokker-Planck equation (|7]) with the initial condition 



U(y,t \yo,t Q ) = 5(y - y ), 



and can be expressed as |30 



(90) 



n(y,t\yo,t c 



-v(y)+v(y ) 



J] $ k {y)$ k {y )e 

k=0 



-A fc (t-to) 



Here, 



v(y) 



4tt 2 



V(y), 



and $>k{y) form a complete ortho normal set of eigenfunctions of the equation 



d 2 



dy< 



+ v'(yf-v"(y) 



My) 



87T 2 A fc 

H 3 



My)- 



(91) 



(92) 



(93) 



The eigenvalues A^ are non-negative, with the smallest eigenvalue Ao = corresponding to 

$ (y) = N-V 2 e- V< »\ (94) 

with 
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N = / e- My) dy. (95) 

J — oo 

Note that $o(y) is related to the "equilibrium" distribution function (|8|), 

P (y) = Sg(y). (96) 

we shall assume that the eigenf unctions are ordered so that < Ai < A 2 < ... . For a purely 
quartic potential [|J, 

A x = 0.0889A 1/2 #, A 2 = 0.289A 1/2 #. (97) 

Let us now consider two co-moving positions x x and x 2 separated by a distance ax = 
o|xi — x 2 | ^> at some moment t = 0. These positions were within each other's horizons 
prior to the time 

t a = -H~ l \n(Hax). (98) 

For t < t a , the values of y at x x and x 2 are essentially the same. Hence, the probability for 
y to take a value y\ at x x and a value y 2 at x 2 at time t = can be expressed as [|J 

P 2 [yi(xi),y 2 (x 2 )] = lu( yi ,0\y a ,t a )Il(y 2 ,0\y a ,t a )Po(y a )dy a , (99) 

where -Po(y) is the equilibrium distribution (j9~6D . The equal-time two-point correlation func- 
tion for y n is given by 

roc 

(^( Xl ,0)^(x 2 ,0)) = (y»y») = / dy 1 ^ 2 y^ 2 "P 2 [y 1 (x 1 ),y 2 (x 2 )] 



= Ar_1 E^(n) fc (^)~ 2Afc/H ; (100) 

fc=0 

where 

A {n)k = / y n e-<y^ k (y)dy (ioi) 



and we have used Eqs. (|9l|) , (0), (0) and the orthonormality of the functions The 
mode functions $fc(y) are even functions of y for /c even and odd functions of y for k odd 
($&(?/) has k nodes), and it is clear from Eq. (|101|) that At n \ k are non-zero only when n and 
k are both even or both odd. 

The first term in the series (|100|) is 

N- l A\ n)0 = (y n )\ (102) 

and thus the reduced two-point function, 

c n {x) = (jMS) - (y n ) 2 , (103) 

is given by the same series starting with k — 1. The asymptotic behavior of c n (x) at large 
x is determined by the first term in that series. For odd values of n, 
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c n (x) « N- x A\ n)1 {Hax)-^' H . (104) 

For even values of n, A( n )i = and 

c n {x)^N- l A\ n)2 {Hax)- 2 ^ H . (105) 

Eqs. (|TT|), (|T2|) of Section III follow immediately from these relations. 
Quite similarly, the three-point function can be expressed as 

(ViVht) = J dyidy 2 dyzylyly -| P 3 [y y 2 (x 2 ),y 3 (x 3 )]. (106) 

Assuming that xi, X2, X3 form an equilateral triangle of side ax, the probability distribution 
P 3 can be written as 

■P3 [3/1, 1/2,2/3] = / R(yi,o\ya,ta)R(y2,o\y a ,t a )n(y 3 ,o\y a ,t a )P (y a )dy a , (107) 



Combined with equations ( |106| ),(|9T|),(p8|), this gives 



00 „ 
(yUyi) = N- 1 £ A (4)fc A (4) ^ (4)m (Faz)-( Al+A2+A3 )^y$ fe (y)^(2/)^(2/)^^^. (108) 



fc,/,m=0 



Once again, it is easily understood that the reduced three-point function (2D) is given by 



the same series with the summation starting at k,l,m — 2. The leading term at large x is 

Ux) « N~ l A\ A)2 C{HaxY^ H , (109) 

where 

/oo 
^( y ) e «(»)dy. (110) 
-00 

^From (105) and ( |109| ) we find that at large x 

Ux)/& 2 {x)^N x l 2 C, (111) 

independent of x. 



23 



REFERENCES 



[1] D. J. Eisenstein, W. Hu, and M. Tegmark, Astrophys. J. in press (|astro-ph/9807130|) , 

and references therein. 
[2] A. Vilenkin, Phys. Rev. D27, 2848 (1983). 

[3] A. A. Starobinsky, in Field Theory, Quantum Gravity and Strings, ed. by H.J. de Vega 

and N. Sanchez (Springer, Heidelberg, 1986). 
[4] A. A. Starobinsky and J. Yokoyama, Phys. Rev. D. 50, 6357 (1994). 
[5] L. Kofman and A.D. Linde, Nucl. Phys. B282, 555 (1987). 
[6] P. J. E. Peebles and A. Vilenkin, Phys. Rev. D59, 063505 (1999). 
[7] G. Felder, L. Kofman and A. D. Linde, |hep-ph/9903350 . 



[8] The one-loop effective potential for the model ([!]) in de Sitter space was calculated 
in Ref. ||. The condition (|J) is obtained by requiring that the one-loop correction is 
negligible for y 2 < (y 2 ) ~ \~ X I 2 H 2 , where the expectation value (y 2 ) is taken in the de 
Sitter- invariant state (see section III). 
[9] A. Vilenkin, Nucl. Phys. B226, 504 (1983). 
[10] E. Gaztanaga and P. Fosalba, Monthly Notices Royal Astronomical Society 301, 
534 (1998); M. White, |astro-ph/9811227| ; J. A. Frieman and E. Gaztanaga, |astro- 
| ph/9903423 . 

[11] A. D. Linde, Particle Physics and Inflationary Cosmology (Harwood Academic, Chur, 
1990). 

[12] L. H. Ford, Phys. Rev. D 35, 2955 (1987). 

[13] P. B. Greene, L. Kofman, A. Linde, and A. A. Starobinsky, Phys. Rev. D 56, 6175 
(1997). 

[14] Ellis, J. et.al, Nucl.Phys. B373, 399 (1992). 

[15] L. Randall and S. Thomas, Nucl.Phys. B449, 229 (1995). 

[16] D. H. Lyth and E. D. Stewart, Phys. Rev. Lett. 75, 201 (1995). 

[17] G. Dvali, Phys. Lett. B355, 78 (1995). 

[18] M. Dine, L. Randall and S. Thomas, Phys. Rev. Lett. 75, 398 (1995). 

[19] T. Damour and A. Vilenkin, Phys. Rev. D53, 2981 (1996) 

[20] M. Dine, W. Fischler and D. Nemeschansky, Phys. Lett. B136, 169 (1984). 

[21] A.D. Linde, Phys. Rev. D53, 4129 (1996). 

[22] G. Dvali and Z. Kakushadze, Phys. Lett. B417, 50 (1998). 

[23] I. Antoniadis, S. Dimopoulos and G. Dvali, Nucl. Phys. B516, 70 (1998). 

[24] L.P Grishchuk and A.D. Popova, Zh. Eksp. Teor. Fiz. 77, 1665 (1979) [Sov. Phys. 

JETP 50, 835 (1979)]. 
[25] A. Linde and V. Mukhanov, Phys. Rev. D56, 535 (1997). 

[26] P. J. E. Peebles, Astrophys. J. 510, 523 (1999); P. J. E. Peebles, Astrophys. J. 510, 523 
(1999). 

[27] K. Coble et al, preprint ( |astro-ph/9902T95|) . 

[28] S. Yu. Khlebnikov and I. I. Tkachev, Phys. Rev. Lett. 77, 219 (1996). 

[29] N. A. Bahcall and X. Fan, Proc. N. A. S. USA 95, 5956 (1998); P. J. E. Peebles, Phil. 

Trans. R. Soc. London. A 357, 21 (1999). 
[30] H. Risken, The Fokker-Planck Equation, Springer- Verlag, Berlin, 1989. 



24 



FIGURES 




log,„ (Hax)^ 

FIG. 1. Numerical realizations of the reduced two-point correlation functions of powers n of 
the field y (Eq. |Tl| ), for n = 1, 2, 3, and 4 from top to bottom at the left-hand side of the 
figure. The dotted lines have the intercept at zero separation and the slope at large separation 
of the Fokker-Planck approximation (Eqs. |l0| to |l2|). The curves are averages across numerical 
realizations of the random process for A = 0.1 (dashed lines) and A = 1 x 10 -4 (solid lines). The 
correlation functions and the separation are scaled to remove the dependence on A and H (Eq. |i~7| ). 




FIG. 2. The reduced three-point correlation function £3 for equilateral triangles scaled by the 
two-point function £ evaluated at the length of the side of the triangle for A = 1 x 1CP 4 (solid line) 
and A = 0.1 (dashed line). In the Fokker-Planck approximation the ratio plotted on the ordinate 
is constant on scales much larger than the coherence length aR c . 
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log 10 x/x x 

FIG. 3. Numerical realizations of dark mass two-point correlation functions at the end of 
inflation for the de Sitter equilibrium state (dashed line) and the rolling Hubble parameter in 
Eqs (|45| ) and (^). The self-coupling parameter in the dark mass potential is A = 0.1 
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